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Abstract 

An individual-based model (IBM) of a spatiotemporal terrestrial ecological popu- 
lation is proposed. This model is spatially explicit and features the position of each in- 
dividual together with another characteristic, such as the size of the individual, which 
evolves according to a given stochastic model. The population is locally regulated 
through an explicit competition kernel. The IBM is represented as a measure-valued 
branching/diffusing stochastic process. The approach allows (i) to describe the as- 
sociated Monte Carlo simulation and (m) to analyze the limit process under large 
initial population size asymptotic. The limit macroscopic model is a deterministic 
integro-differential equation. 
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1 Introduction 

Our aim is to present a spatially explicit individual-based model of a spatiotemporal ter- 
restrial ecological system. We consider a family of individuals whose state includes their 
position and another characteristic such as their size. During the life of an individual, 
its position remains constant and its characteristic evolves according to a given stochastic 
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continuous model. Each individual is subject to random punctual mechanisms: natu- 
ral mortality, mortality due to competition, and reproduction. The individuals interact 
through the competition mechanism. 

So-called individual-based models (IBM) are not new in the field of theoretical ecology 
where they have become more widespread thanks to the development of computers. Since 
the publication of the significant paper by Huston et al [18], many studies have been 
devoted to these models [8j, [17J. IBMs appear mainly as computer simulators based on 
empirical rules. IBMs are modeling tools that require design work, and which cannot 
be reduced to the indiscriminate use of empirical rules. The study of the behavior and 
properties of IBMs still requires mathematical analyses. 

Mathematical representations of such models in continuous time and space was in- 
troduced by Bolker-Pacala j2] and Dieckmann-Law ||9j. These authors also derived a 
truncated moment method that determines the time evolution of the approximated first 
moments. Meleard and co-workers [15], [5]) [1|) [S] proposed a rigorous mathematical setup 
for these models: the dynamics of the ecological system are described as a measure- valued 
branching/diffusing random process. The limit law behavior of this type of microscopic 
Markov process, as certain parameters such as the size of the initial population tend to 
infinity, is relatively well known. The macroscopic limit models can be of different nature: 
deterministic, such as ordinary, integro, or partial differential equations; or probabilistic, 
such as stochastic partial differential equations or superprocesses [E], |13j . 

Within this framework, various kinds of dynamics have been considered: terrestrial 
plants [12], Darwinian evolution [SI HIE]) phytoplankton aggregation [T^, age-structured 
population [23] . 

Meleard and co-workers also described an exact Monte Carlo procedure to simulate 
the associated microscopic stochastic process. These algorithmic aspects which are very 
important in practice, deserve specific attention. 

The objective is threefold. Firstly, we set out an IBM mathematical framework for 
ecosystems such as terrestrial plant systems. Secondly, we develop the associated Monte 
Carlo algorithm. Thirdly, we derive the macroscopic behavior of the IBM. 

We focus on terrestrial plant ecosystem dynamics models that are spatially explicit 
with an explicit representation of the competition interactions between individuals. This 
is one of the most active areas of computational ecology [16], [21], [1]. 

In the authors present a model for terrestrial plant ecosystem dynamics where 
they consider only the position of the plant individuals. In our work we extend this model 
further to include a continuously evolving characteristic such as the size of the individuals. 

We describe the Bolker-Pacala-Dieckman-Law model in Section 2. The Monte Carlo 
simulation method is presented in Section 3. An example of a terrestrial plant ecosystem 
is presented in Section 4. The Markov representation of the IBM is described in Section 
5. Finally the large population limit is analyzed in Section 6. The generic mathematical 
approach used to prove the weak convergence of measure-valued process in Section 6 is 
detailed in the Appendix. 
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2 The model 



We consider a family of individuals that live in a set of the form: 

where P is a measurable connected subset of M". The state x = {p,r) ^ X of an individual 
represents its position p in the physical space P and an associated characteristic r that 
could represent its size or its maturation age. 

It is convenient to represent an individual at point x £ X as the Dirac measure 5x', 
hence the population at time t will be: 

iyt{dx) = ^5,.{dx) (1) 

i=l 

where the sum is over all individuals alive at time t and Nt is the size of the population 
at time t. 

Let A4f{X) be the set of finite positive measures on X, and Ai be the set of finite 
point measure on X, i.e. 

M = {Ef=i^x^;Nen*,xiex}. 

The duality operator between the measures and the functions will be denoted: 

(^,/)= / f{x)Kdx) 
Jx 

(note that {Z)^ will denote the quadratic variation of a process Yt). Hence, (ft)t>o is a 
-valued process and, according to ([l]): 

is the size of the population at time t. By abuse of notation, "x G ff" will specify that an 
individual in the state x belongs to the family vt at time t. 

Remark 2.1 (numbering convention) In practice, starting from a family labelled from 
1 to N, we use the following numbering of particles: (i) in case of birth, the new individual 
will be labelled N+1; (ii) in case of death of the individuali, the first (i — 1)*^ labels remains 
unchanged and the last [N — i)*'* labels are left shifted (i.e. j — )• j — 1 ). As pointed out in 
fTS^ and J^, this numbering convention has no influence on the law of the process we will 
describe, it affects only the trajectorial realizations of the process. 

Considering the state u = YliLi of the family at a given time, an individual in state 
X G v will be subject to 3 types of punctual events occurring at specific rates: 

Intrinsic death: This individual disappears at a rate X'^{x) which may de- 
pend on its state. This death is called "intrinsic" as it does not depend 
on the state of all the population v. It represents the "natural death" as 
opposed to "competition death" . 
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Competition death: This individual disappears at a rate X''{x, u) whicli may 
depend on its state x and on the state v of the population. We suppose 
that \''{x, v) is of the form: 

>^'{x,v) = ^u{x,y)= j u{x,y)v{dy). (2) 

The competition kernel y) is the contribution of an individual located 
at y to the competition affecting an individual located at x. 

Birth and dispersal: This individual gives birth to a new individual at a 
rate A''(a;) which may depend on its state. The state y £ X oi the new 



individual will be determined by a given dispersal kernel (see Remark 2.2 
later) . 

Between discrete events of birth or death, the size of the population remains unchanged 
as well as the position of the individuals, and the population state {x\)i<i<_N is subject to 
a continuous mechanisms: 

Displacement: Over time, the characteristic component of each individual i 
evolves in the state space in interaction with the evolution of all other 
individuals according to the following system of stochastic differential 
equations (SDE): 

d(:i) = (5(4,^.))di+(.(.U))dBj (3) 

where (B^)j>o are independent standard Brownian motions. To simplify 
the notation. Equation ([s]) will be represented as: 

dx\ = g{xiut)dt + a{xivt)d-Q\. (4) 

Let 

a{x, v) = a{x^ v) a*{x, v) . 
We define the associated flow operator: 

vt = ¥{t,s;vs) 

defined for all s < t between two successive punctual events (i.e. between 
two successive jump of the population size). Note that the initial condi- 
tion distribution law in SDE ([3| or Q is handled by the the dispersal 



kernel (see Remark 2.2 later). 



We suppose that these four mechanisms and the dispersal mechanism are mutually inde- 
pendent. 
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Remark 2.2 (dispersal kernel) An individual in state x = {px^fx) will give birth to a 
new individual in state y = {py,ry) = x + y' = {px+Py' ^fx + T^y') ("y" o,nd "y' " will denote 
respectively the absolute state and the relative state of the new individual). The state of 
the new individual is given by a relative kernel. 

In the present application it will be convenient to consider a "mixed relative/absolute" 
formulation: an individual in state x = {px,rx) will give birth to a new individual in state 
(px -\-Pz,Tz) according to a kernel D{x,(iz) with z = (pz^fz). We suppose that this kernel 
admits a density: 

D{x,dz) = D{x,z)dz . (5) 

This formulation is natural: the position px + Pz of the new individual will be relative to 
Px and its characteristic r^ will be absolute. For phenotypic trait dynamics l^jj, it is more 
natural to consider a relative mutation kernel for the r-component. 

3 Monte Carlo simulation 

We now describe the dynamic of the process starting from an initial population state v. 
Independently of each other, an individual with state x in the population v has three 
independent exponential clocks that control the occurrence of the events: (i) a birth clock 
with rate A''(2;), (ii) an intrinsic death clock with rate A'*(x), (Hi) a competition death 
clock with rate X''{x, i^). 

(i) When the birth clock rings, the individual x gives birth to a new individual with a 
state z G X determined by the dispersal kernel D(x, dz): 

(a) When the intrinsic or competition death clock rings, the individual x is removed 
from the population. 

v ^ u — 6x ■ 

Between any birth or death event, the state of all the population evolves according to 
Q which corresponds to a system of A'^ interacting SDE's (A^ = {u, 1) is the size of the 
population). 

Considering individual clocks is cumbersome, a more efficient Monte Carlo procedure 
will rely on the existence of a global clock that dominates all punctual phenomena (birth, 
natural death, competition death). That existence holds true when all the different local 
clocks are uniformly bounded: then, given an individual chosen at random in the popu- 
lation, the type of punctual phenomenon to be considered is determined by a sampling 
technique, and it is decided whether the chosen phenomenon is actually applied or not 
by an acceptance/rejection sampling technique. The existence of a uniform bound avoids 
explosion phenomena, i.e. accumulation of infinitely many events at a given time. 
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Hypotheses 3.1 We suppose that there exists positive real numbers A^^^, A^^^ and u 
such that: 



X\x) < Al, , X'ix) < Al, , u(x, y) < n_ . 

Hence: 

A'(x,z^) < U,^aA^A) ■ 
We also suppose that the dispersal kernel ^ satisfies: 

D{x,z)<kD{z) (6) 

where D is a probability density function. □ 

Let To = 0, and suppose Tk-i and ^T^.i given. We describe now how to simulate ut^ 
starting from i^Tk-i- order to determine the instant where the next event could take 
place, we should bound the different rates uniformly in space but also in time. This is 



possible thanks to the Hypothesis 3.1 From the instant Tf^-i to the instant of the next 
event, i.e. along the time interval [Tk-i,Tk[, the population size is = {i^Tk-i^ !)• 

At the scale of the population, the maximum rate of events (birth, natural death, 
death by competition) is bounded by: 

7 = 7*^ + + 7<= with { 



I max 

7^' = A<^ N (7) 

/ max ^ / 

Y = u N"^ 



where 



7*^ is an upper bound of the birth rate at the population scale (i.e. the birth 
rate at the population scale if all individuals have the same birth rate 

\b \ 

max / * 

7^^ is an upper bound of the natural death rate at the population scale (i.e. 
the natural death rate at the population scale if all individuals have the 
same natural death rate AjJ^^^). 

7^^ is an upper bound of the rate of death by competition at the population 
scale. 

An acceptance/rejection method will permit us to correct the fact that these three terms 
are upper bounds for the actual rates. 

Iteration I'Tk-i ~^ '^T^- 

(i) Let N = {uTk-i^ 1) be the population size, 
(a) Computation of the global rate 7 == + 7*" + 7^" with 
(Hi) Simulation of the next event instant: 

Tk = Tk-i + Sk , with Sk ~ Expi'j) . 
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(iv) Computation of the system evolution between the two instants: 

Uj,- = ¥{Tk, Tfc.i, i/Tfe.J . 

In practice, the system is simulated with an Euler discretization scheme. 

(v) Chose X at random uniformly in Vrp- ; chose at random the nature of the 

k 

next event according the probability values (7''/7, 7'^/7, tVt)- 
• birth: choose z' according to the law D{z)dz and let 

X>'{x) D{x,z') 



{P.+P..,r^,) w^ih probability ^//fff.^)^ 



Vrp- + (5, 

Vrr- with probability 1 — ''^i^'* '^fc'^,? 

• natural death: 

= 



Urp- — 5x with probability \d^^ 
Vj,- with probability 1 



• competition death: chose y at random uniformly in Urp- and 

k 

let 



Vrr- — 6x with probability "^^'^^ 



Vrp- with probability 1 — "^^'^^ 

k 



□ 



The algorithm is detailed in Algorithm [T} 

We will present numerical tests in a future companion article [3] where we will also 
propose other versions of this algorithm. 

4 A forest dynamic model 

Competition model: zone of influence approach 

Here we consider a population of trees. The state of each individual tree i is: 

X* = {p\r') e [0,Lf X [r'-^r"^^"] C x M+ 

where is the position of the tree in a parcel V = [0, L]^ and r* is the radius of its zone 
of influence (cf. Figure This zone of influence is the disk centered in with radius 
which symbolizes the portion of the ground that the individual needs to ensure its growth. 
For X = {p,r) £ I' let: 

Vx = disk of center p and radius r. 
The local interaction kernel u{x, y) introduced in ([2]) is of the form: 

r Arca(P.nP,) -r / 

U{x, y) = I Area(C.) " ^ y , 

otherwise. 
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1 To ^ 0, /c 0, initialization of 

2 while Tfc < Tn^ax do 



3 
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k k + l 

6 5 ~ Exp(7) , Tfc ^ Tfc_i + S" 

7 Urp- ^ ¥{Tk,Tk-i,h'Tk-i) {evolution (Euler scheme)} 

8 u^U[0,l] 

9 choose X at random uniformly in Vj,- 

10 if n G [0, 7''/7] then 

11 v^U[0,l] 

12 z' ~ D{z) dz 

13 if < (X'^ix) D{x, z'))/{Xl^^K D{z')) then 

14 z^Tfe ^ J^T,_ + <5(p,+p,/,r,,) {birth} 

15 end if 

16 else if u G [7''/7, (7'' + 7'^)/7] then 

17 vr^UlO,!] 

18 if u < A''(x)/A'^ then 

— V / / max 

{natural death} 

20 end if 

21 else 

22 choose y at random uniformly in Vj,- 

23 i;~ZY[0, 1] 

24 if w < n(x,y)/Mmax then 

25 VTf, z^T _ — {competition death} 

26 end if 



end if 



28 end while 



Algorithm 1: Simulation algorithm, first version (see description in Section^. 
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Figure 1: Left: Each individual tree i is associated with a state = {p^,r^) where is 
its position m D C and r* G [7-^'°^ T-m^xj radius of its zone of influence. 

This area represents the portion of ground which the individual needs to ensure its growth. 
Right: The more the zone of influence T>x of an individual in state x intersects with other 
individuals ZOI, the higher the strength of competition is, and the greater the associated 
death rate is. The ZOI appears in the definition of the competition kernel u{x, y) in ^ 
and in the growth model (10). 



The surface area Area(Pa;) of the zone of influence associated with an individual in state 
X in z/ represents the amount of resources (e.g. sunlight, water, nutrients...) needed for 
growth; u{x, y) is the strength of competition experienced by the individual in state x 
from an individual in state y in the population u. 

The computation of the areas in ^ could be cumbersome for large population sizes. 
However, there are many alternatives pj. 



Birth and dispersal model 

Birth occurs at a rate X^{x) defined by: 

A^(x) = A>, r) = AL. ^ l{.>.n < • 

This birth rate can be understood as a fertility model: only individual with ZOI radius 
greater than r^ can give birth; and the greater this radius is, the greater the rate is. 
For the dispersal kernel D{x, z), we can consider two possibilities: 

(i) "Parcel in forest model" - Here T> is the torus [0, L]^, and we consider an homoge- 
neous kernel: 

D{x,z) = D{z)=Dx{pz)D2{r,). 

In this example we consider a slightly different case where the initial condition on 
does not admit a density: For example: 

D{x, dz) = D{x, dpz X drz) = M{0, cj^ I; dp^) x 5rmin{drz) 
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where is a "Gaussian law on the torus P" 

The condition ([g]) is fulfilled but in fact the present setup is even simpler. Indeed, 
the part (v) of the algorithm proposed at the end of sectionals now: 



• birth: 



^T' + ^(p^+p,^,r,,) with probability ^ 



i^j,- with probability 1 



where z is simulated according to the law Di{pz) D2{rz)- 

This setup is periodic and could illustrate the case of a squared parcel imbedded in a 
forest: the descendants of the individual of the parcel may integrate the contiguous 
parcels and the individuals in the contiguous parcels may integrate the parcel under 
consideration. 

(ii) ^'Island model" - Here P = [0, L]^ is a bounded squared parcel and we suppose that 
when an individual is closed to the border of P a portion of its offspring is lost in the 
water and do not integrate the population. Hence, the individuals near the border 
have a lower effective fertility rate compared to the individuals closed to the center 
of the parcel. 

For simplicity we suppose that the birth fertility rate is constant for individuals 
at distance greater than ro from the border, for these individuals suppose that the 
dispersal kernel is homogeneous: Di{px,Pz) = Di{pz) (with the support of Di{pz) 
is included in a disk of center and radius rg), suppose that all these individuals 
have a constant fertility rate A''. For individual located in px at most ro from the 
border, then 

Dl{px,Pz) =CxX l[0,L]2(Px +^2) X Di{pz) 

with Cx = [/[Q 2,]2 Di{py — Px) dpy\~^ . Here we suppose that the fertility rate associ- 
ated with the later points is non-constant and equal to \°{x) = X° /Cx, i-e. the more 
the kernel Di intersects the border, the lower the birth rate is. 

In the previous model we compensated the birth rate for individual close to the border. 
Without this mechanism we get a model where there is an accumulation of new individuals 
on the border of the parcel. In a way, the "parcel in forest" model is the more realistic 
and therefore, this is the one we used for simulation purposes. 



Growth model 

Suppose that the radius of the ZOI is solution of the deterministic equation: 

rl=9(xlvt) (9) 

^That is M{Q,a'^ I;B) P(X G Up.qezB + {pL,gL)) for any Borelian set B of [0,Lp where X ~ 
AA(0,a2j). 
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while the position remains unchanged. This last equation is coupled with = so that 
we get a system of ODE's for xj. For ([o]), we consider a model of the form: 

ri = i,{xi,ut)R{ri), ro = r^'" (10a) 

where R{rl) corresponds to a standard growth model such as the Richards model [7]: 

R^rr=ai^.Yh^j\[{^r-'-i] (lob) 

with 7^ 1 and ip{x\,vt) is an expression between and 1. The case tjj{x\,vt) = 1 
corresponds to the best condition for growth: in this case its growth is described by the 
Richards model r\ = R{rl). The smaller ip{x\,i't) is, the more the growth conditions 
of the individual i are degraded. Small i{j{xl,iy't) corresponds to the situation where the 
individual i is surrounded by many other individuals. We may think of: 

iPixl, ut) = [l - Cg X^{xi, ut)] ^ G [0, 1] (10c) 

or, more generally, il){x\,ut) = i^'' vt)) where ^ : M+ i— )■ [0,1] is any continuous, 
decreasing function such that ^'(0) = 1. 

5 Markov representation of the process {yt)t>{) 

5.1 Identification of the infinitesimal generator 

We introduce the following set T) of test functions ^ : Mf{'^) i— M of the form: 

^u) = F{{vJ)) 

for any function f : X ^ and F : M i— ?• M twice continuously differentiable, bounded 
with bounded derivatives. 

At the end of this we will present a particular case of such test functions ^. 

Lemma 5.1 (expression for ^{yt)) For any <I> = {F,f) G V, ^{vt) = F{{''^t,f)) satis- 
fies: 

^{ut) = ^{uo) + [ C^ii^s) ds + M$,t(i/) 
Jo 

= $(^o)+ Yl { f C^Hi^s)ds + Ml,iu)} 

A=d,b,c,g ''^ 

where 
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(i) C = + + C + are the infinitesimal generator defined by: 

jC'^u) = f X%x) - 5^) - ^{u}] u{dx) . (11a) 

/(dx) (lib) 

/ u{x, y) u{dy)\ [^{v - 5^) - ^{u)] u{dx) (11c) 

C-^ii.) = F'{{u, /)) X {u, gf{.,u)) + lF"{{u, /)) X {u, |V/ • a{; z.)|2) (Ud) 

corresponding respectively to natural death, birth, death by competition and growth 
respectively, and with: 

gf{x, u) = V/(x) • g{x, + l £^ aeeix, v) 

{ii) the martingale term M$^f(z/) = ^(y) + ^(y) + ^(y) + M|, ^(z/) is the sum 
of Jour martingales defined by: 





JX 




[ A 




JX 




def 1 




JX - 







X [$(i/,- -5^i_) -$(z/,-)]N''(ds,di,d^), (12a) 

s 

^lA") = Jo J^,J.^Jq ^(^<^.-) l{e<(A^^_)D(a:;_,^))/(A^,,«D(z))} 

X [$(i/,- + 5(p^,_+p„r.)) - Hi^s-)] ^\ds,di,dz,d9) , 

(12b) 

X [^{u,- -d,i ) - $(1/,-)] N%ds, di, dj, d0) , (12c) 

s 

M|,,(i^) = /* F'iius, /)) X V/(a;|) • a(x|, z.,) dB^ , (12d) 

JO 

and 

a. N"* is a Poisson random measure on [0, oo) x N* x [0, 1] of intensity measure: 

n^{ds,di,de)'^\i^dsdide 

( di is the counting measure on W , ds and dO are the Lebesgues measures on 
[0,oo) and [0, 1]), and N'^ = N"* — n'* is the compensated measure. 

b. N*" is the Poisson random measure on [0, oo) x N* x ^ x [0, 1] of intensity 
measure: 

n''(ds, di, dz, dO) = Xl^^ kD{z) ds di dz dO , 
and N'' = N** — n'' is the compensated measure. 
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c. N'' is the Poisson random measure on [0, oo) x N* x N* x [0, 1] of intensity 
measure: 

n''{ds, di, dj, d9) = u,^^^ ds di dj dO , 
and N"" = N"" — n'' is the compensated measure. 

d. (Bj) 

t>o, i>l, is a family of mutually independent standard Brownian motions 
in R. 

These four stochastic processes are mutually independent. 

Proof As the four basic mechanisms (natural death, birth, competition death, growth) 
are independent, we can write: 

^ut) = ^uo) + Tt + T^t+Tl + TI (13) 

where T^,T^, Tl, r| are the terms associated with the natural death, the birth, the death 
due to competition, and the growth. We consider the four terms successively. 

Death component 

The Monte Carlo procedure described in Section [3] implies that: 

I LI ^^'-"^^"^ l(e<Ad(.»^_)/Ad^j [H'^s- - S^iJ- H^s-)] N<^(ds,di,de) 
where N'^ is a Poisson random measure described in [ii)-a. Here we use the labeling con- 



vention given at Remark 2.1 , is the i particle of Vg- . By introducing the compensated 



measure N"* = N** — n'^ and the martingale M|, ^{v) defined by (12a), we get: 

r? = j^^ \^<N^^) l(e<Acl^.)/Af_) [H^s- - K^J- H^s-)] AL.dsdid0 + M.%^^{l 
= f I \^<N^.)>^'K~)M^s~-K^ )-^U,-)]dsdi + m%^,{u) 

Jo Jm* 

= / )-Hu,-)]ds + M%/u) 

= [ [ X'ixJMu,- -6^)-^iy,^)]iy,^{dx)ds + M%/u) 
Jo Jx 

Birth component V\ 

The Monte Carlo procedure described in Section [3] implies that: 



/q LIxIq hi<N,^)'^{e<(\^{xl_)D{xl_,z))/(\^^^^^D{z))} 



X [^.(z.,_ ^ +p^,,^)) -$(i/,-)]N'^(ds,di,dz,de) 



13 



where N'* is the Poisson random measure described in (ii)-b. By intro ducin g the compen- 
sated measure N*^ = N*" — n'' and the martingale M|, ^(i/) defined by (12b), we get: 



JN* JX JO 



l{i<7V^_ ) l{e<(Ab(x'- _ ) Dixi_ ,z))/iX^,^ K D{z))} 



X [^V,- + 5. 



j) -*(■'.- )1 AL. K d» di dO + Mi/f) 



))-cI>(z.,-)]Z)(x:_,z)dz} ds + M|,,, 



Competition component 

The Monte Carlo procedure described in Section [3] implies that: 



l(i<7v^_) l(i<7v,-) l(e<«{x'_y _)/«„,ax) 

\ s s ' ' ' 



JN* JW JO 



X -(5,. )-$(i/,-)]N^(ds,di,dj,de) 



where is the Poisson random measure described in {ii)-c. By intro duci ng the compen- 
sated measure = N"" — n'' and the martingale M|, ^(zv) defined by (12c), we get: 



JN* JN* Jo 



l(i<JV^-) l(i<^,-) l(e<«(a-_y_)/«^ax) 



JN* JN 



s 

X - 6^. ) - )] ds di dj + Ml , (i/) 



.N- N 

ft a 



Growth component Vj. 

Consider now any instant between two instants of jumps (i.e. between two discrete events), 
for simplicity's sake we consider a time interval of t G [0, to] between and the first instant 
of jump. For t E [0, toli the evolution of ut is modeled by the SDE (|4]), that is: 

dxi = g{xi, vt) dt + a{x\, vt) dB\ . 

For all ^ = (F, /) G V, from the Ito formula: 

/(^9 = /(4)+ f Qf{xl,ys)ds+ fvf{xl)-a{xi,Us)dBi 
Jo JO 
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with the summation convention for repeated indices and where a = a a* . Hence: 



where 



{vtJ) = {v^J)+ / {us,Qf{-,iys))<is + Yt 
Jo 



ft Ns 
■^0 i=l 



This last expression is a martingale with quadratic variation: 

ft ft 

{Y)t= j;|V/(x:)-a(xl,z.,)|'ds= / (z.„|V/-a(-,i.,)|2)ds 



We get: 



F[{ut,f)) = F{{u^J))+ / F\{us.f))>^Mf{;ys))ds 

Jo 

Jo 

t 



+ 



F"{{usJ)) X (z.„|V/-a(.,i.,)|2)ds 







The infinitesimal generator associated with the growth phenomenon is obtained by taking 
the expectation of the previous expression. Finally 



Jo 



where M|,^(zv) is the martingale defined in (12d). 



Summing up these results in Equation (13) ends the proof of the lemma 



□ 



As a corollary we now compute the predictable quadratic variation of the martingale 
processes Mj^(z^): 

Corollary 5.2 The predictable quadratic variation associated with martingale processes 



M;'($), M;'($), M;^($), Ms($) defined in Equations (7^ are: 

{m%^{v))t= f I \\x)[^u,--5,)-^v,-)fu,-{dx)ds, 
Jo Jx 

t 



(14a) 



{m%^iu))t= [ [ X\x) [ Mu,^+6(^p^+p^^r.^)-^u,-)fD{x,z)dzu,-{dx)ds, 
Jo Jx Jx 

(14b) 



m>M))t 



[\[ u{x,y)u,-{dy) mu,- -5^^)-^{iy,^)]'^u,-{dx)ds, (14c) 
Jx '-Jx 

F'i{u,-,f))\'\Vfix)-a{x,u,-)\\,-idx)ds. (14d) 

Jx 
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and by independence of the processes N"*, N'', and B*, z > 1, we get: 

Proof Consider a martingale process: 

Mt= [ [ K(s,x)N(ds,dx) 
Jo Je 

where N is the compensated measure associated with a Poisson random measure N on 
[0,cx)) X of intensity measure n (i.e. N = N — n), and k{s,x) is a predictable process 
such that E ^^(s, x) n{ds, dx) < oo. Then: 



{M)t = [ [ K^{s,x)n{ds,dx). 
Jo Je 



(see [191 Section II-3]). We apply this result for example for the first term, from (12a 

ft 



JO JN* JO " 

X [$(z^,- -6^. )-<^>{i^,-)f n^(ds,di,d0) 



= I J^,l hi<N,-) l(e<Ad(x^_)/AjLax) - ^^iJ- '^(''s-)? Al^d^dids, 

integrating in 9 leads to 



(M|„(^))*=/V l(,<iv,_)^5^[$(i^.--5.. )-$(i.,-)]2AL.did. 

Jn* ^ - " ' ''max 

V \^<N^)>^^'{xl-)ml^s--S^' )-Hu,-)]^dids 

JN* 

X-'ix) [^u,^-6,) - ^,^,-)]^u,^{dx)ds 

JX 



The same approach could be applied to the next two expressions, the last assertion (14d) 
is due to the properties of the Brownian motion. □ 

We now consider a particular case of test functions $: 

<^(^) = ^,/) 

(i.e. with F =id.) for any function / : i— )• R twice continuously differentiable, bounded 
with bounded derivatives. 
We define: 
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for A= d,b,c,g. Hence: 

{i^u f) = {i^o, f) + f ^fi^s) ds + mf^t{u) 
Jo 



Wo,f)+ Yl { A"/(^s)ds + mf,,(z.)} (15) 



A=d,b,c,g 

where 

(i) i = i"^ + + + i'^ are the infinitesimal generators defined by: 



£VM = - / A^(x)/(x)Kdx), 

tf{v)= [ X^{x)[f^f{ip,+p,,r,))D{x,z)dz]u{dx) 
Jx 

= - / [fx y) ^(dy)] fix) iy{dx) , 

^VM= / Gf{x,u)v{dx) 



X 

corresponding respectively to natural death, birth, death by competition and growth 
respectively. 

[a) the martingale terms mj((i/) = mjj(i/) + m°j ^{v) + m^^ ^{v) + m^j^{v) are defined 
by: 

= - j^^j^ 1(*<^,-) l(e<Ad(x'^_)/Ad^j /(x*-) N^'(ds,di,d^) , 

= ^ j^JxIo ^^'-^^'^ l(e<Ab(x'^_)/Ab /(a^s- + ^) ^''(ds, di, dz, d0) , 



Jo Jn* Jn* Jo 

X /(4-)N^(ds,di,dj,de), 

where the compensated Poisson random measures N'^, N'', N'' and the Brownian 



motions (Rl)t>o, i > 1, are introduced in Lemma 5.1 The corresponding predictable 
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quadratic variation terms are: 

{nil,{u))t = fix) M<ix) ds , (16a) 

{in}^Xu))t= [ [x^{x)[J;^f{{p^+p„r,))D{x,z)dz]iys{dx)ds, (16b) 



/o Jx 

{uilXu))t= [ [ [f^u{x,y)i^s{dy)]f{x)us{dx)ds, (16c) 
Jo Jx 



{inl{u))t = J JjVf{x)-a{x,Us)\\s{dx)ds. (16d) 
and by independence of the processes N"^, N'', N'' and B*, i => 1, we get: 

5.2 Control of the size of the population 

Lemma 5.3 If for some p > 1, IE((fo, 1)^) < oo, then for any T > 0, 

E sup {i/t, 1)^ < oo . 

0<t<T 

Proof We apply Lemma |5.1| with -F(^) = and f{x) = 1. As the terms corresponding 
to death (natural death, competition) are less than zero, and as = 0, we get: 

iVf < < + / X. X lo ^'-'-'"'-^ ^(^<^'^K-)AS.aJ [(^- + 1)'' - N'(ds, di, dz, de) . 
We introduce the stopping time: 

Tn = mf{t >0;Nt>n} 

then 

X N'^(ds,di,dz,de) 
as (A; + I)*' — k'^ < Cp {1 + kP~^) for some constant Cp, we get: 

l-T/\Tn I- I- l-l 

sup Nf<Ni; + Cp / / / l(,<^_)[l + iVjr^]N''(ds,di,dz,de). 

o<t<TATn Jo Jn* Jx Jo 
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Taking expectation leads to: 



E sup iVf < E(iV^) + CpE [ "Iff l(i<N _) [1 + N^-^] n'(ds, di, dz, dO) 

Q<t<TArn Jo Jn* JX Jo 

< E{N^) + Cp r^m.^TAr^^ [1 + NP-']) ds 

Jo 

< E«) + Cp [ ^lis<TAr^) [1 + N!]) ds 

Jo 

<K{NI) + 2CJ E(iV,Vjds 

^0 

<E{NP) + 2Cp[ E( sup <ArJdt (17) 

JO 0<s<iAT„ 

and by Gronwall's lemma 

E sup Nf < Cp^T (18) 

0<t<TATn 

for all n. We new want to check that r„ — )• oo a.s. Suppose that 7^ 00, then there 
exists To < 00 such that eo = IP(sup„ r„ < Tq) > 0, then: 

E( sup Nj^)>nPeo 

0<i<ToAr„ 



which is a contradiction. Letting n — )■ 00 in (18), by Fatou's lemma we prove the result. 
Note that in Equation (17) the different constants Cp depend only on X]"^^^ and on 



universal constants, but do not depend on the functions u, g, a. □ 

6 Large population limit 

Let k be the initial population size, i.e. k = {uq, 1), and replace u by u^, g by g^ and a 
by (T*^. Let {iyt)o<t<T be the Markov process defined in the previous section with initial 
population size k. We define: 

I. def 1 k 

In this section we study the asymptotic property of the law of the process (/ij )o<t<r on 
the space O{[0,T],J^f{X)) of cddlag functions from [0,r] with values in J^f{^) (Q- 
According to the general approach depicted in the appendix, a classical method to do this 
consists of deducing both convergence and characterization of the limiting process from 
the convergence and properties of (/x^ , /), for all / in a suitable class. 



From (15): 



{u^J) = {vlf) + fi'fiu^) ds + in%{u') (19) 
Jo 



^If not mentioned "BdO, T], Mp{X)) is equipped with the Skorohod topology associated with the weak 
topology on Mf (X) see appendix. 



19 



where: 

i'fii^) = J^{-\'{x) f{x) + X^{x) [ f{{p. + p., r,)) D{x, z) dz] 
- u{dy)]f{x) + g''f{x,u)]u{dx) 

with 

g'fix, u) V/(x) • gHx, u) + i ^££r; 4,{x, v) 

and 



/O JN* JO 



m^,tM = -/ / / l(*<w,-)l(e<Ad(x»_)/Ad_)/(a:l-)N'^(ds,di,d^) 

JO JN* ^0 

l(i<jv^_) l(e<Ab(x^_)/Ab^j f{{Vx\_ +P^,r^))N'^(ds,di,dz,d6') 
l(i<JV,-)lo-<7V^-)l(e<,fc(,^ y N^(ds,di,dj,d0) 

+ (Y.tx^fi.<)-<y\<.ys)dw, 

Jo 



10 JN* JX JO 

r-t r f fl 



JN* JN* Jo 



Dividing (19) by A; leads to: 



Define: 



and 



with 



{„l f) = {„l f)+ [\ ^'fi'^s) ds+l m^ii.') 

= it^i f) + [\ ^'fikf^s) ds + l mJ,,(A;/) (20) 

Zr=lm%ikf^') (21) 



r + r + t + (22a) 



r/(/i) = ^r/(fc/i) = e'fifi) = - \\x) fix) /i(dx) , (22b) 

tfifi) = ^tf{kfi) = tf{i^) = J^^'i^) [f;^f{{p.+Pz,r.))Dix,z)dz] is{dx), (22c) 

^"7(/^) = - / [ ^(^' y) /^(dy)] fix) Kdx) (22d) 

JA- 

^V(/^)= / Gf{x,fi)fi{dx) (22e) 
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AT 



and 



gf{x, n) = Vf{x) ■ g{x, + 5 £^ aw {x, 



Hence (20) reads: 



where 



{^^h f) = {^^o, f) + / mf^s) ds + zf + Rf 



k 



ds. 



From ( 16 ) we get: 

= \[ {a-^I^) P{x) + \^{x) [ + p„ r,)) D{x, z) dz] 

+ |V/(x)-a'=(x,fc^J)|2}^,fc(dx) ds. 

We are now ready to state the main result of this section. 
Theorem 6.1 Suppose that: 

k u (x, y) — ; > u{x, y) , 

fe— >oo 

y v-f-) ^ y\x, , 



L°°(A'xXp.(A')) 



fc— ^-oo 



and t/iai g{x,C} and a{x^^) are hounded and Lipschitz in x uniformly in ^, i.e. 

\9{x,0 -9{x,^')\ + \a{x,0 - a{x,^')\ < C p{fi,fi') 



(23) 



(24) 



(25) 



(26a) 
(26b) 
(26c) 

(27) 



where p is the Prohorov metric on Mf{X) (p generates the topology of weak convergence 
on Mf{X)). Suppose also that: 



h law J. 

/^S > ^0 . 



(28) 



where G Aip^^) (deterministic). 

Then {p.^)k£N converges in law to a deterministic process ^ G C{[0,T],^AF{X)), char- 
acterized by 



(29) 



{^t,f) = {^o,f)+ / m6)ds 



where if is defined in Equations (22). 
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Remark 6.2 Hypothesis (26a) can he understood as a "small competition /la rge popu- 
lation" asymptotic. In the example of Section^ the term is given by (10) and the 
convergence (26b) can be deduced from (26a) as: 



/(x, kfi) = ^{Jk u^(x, y) fi{dy)) R{r) > ^ (/ u^{x, y) fi{dy)) R{r) 

k—>-co 

Note also that in this example a^{x, k jj) = 0. 



The end of the section is devoted to the proof of Theorem 6.1 



Uniqueness 



Lemma 6.3 Equation (29) admits a unique solution in C([0,T], 



Proof First note that, as for the proof of Lemma 5.3 

rt 



(6,l)<(Co,l)+ / / A'=(x)6(dx)ds<(Co,l) + AL. / (6,l)d^ 
Jq Jx Jo 

so that, from Gronwall's lemma, sup^gjg 2^](^t, 1) < ('^o, 1) 6xp(A['j^3^ T) =: kt- 



Consider two solutions and of (29) with the same initial condition we prove 



that W^t — CiIItv = for all t G [0, T] where: 

||/i||Tv= sup \lJ,{B)\= sup |(/x,/)|= sup 

BeB{X) / Borel function f&Cg°(X) 

II/II-<1 ll/l|oo<l 

where Cq^{X) is the set of smooth functions with compact support. For any / G Cq^{X) 
with ll/lloo < Ij we have: 



I (e* - L f)\< f^{\j^ X'ix) fix) Mdx) - 6(dx)] 

( fiip. + Pz, r.)) D{x, z) dz) [^.(dx) - 6(dx)] 



+ 
+ 
+ 
+ 
+ 



X 



Jx "(^' y) ^s{dy)) f{x) [6(dx) - 6(dx)] 



u(x, y) [6(dy) - U^y)]) f{x) 6(dx) 

X 

(V/(x) • + \ 0^7 -au'{x,is)) [6(dx) - 6(dx)] 

(V/(x) • [g{x,is)-g{xA3)] + \ gjoli [a«/(x,Cs) - aee' 



X 



X 



(x,6)])6(dx)|}ds 
(30) 
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The first five terms of the RHS of ( 30 ) are controlled the same way, for example the third 
term is: 



<i 



where u{x,y) < limax- For the last term of the RHS of (30) we use the fact that the 
coefficients g{x,^) and a{x,(,) are bounded Lipschitz in x uniformly in ^. We get: 

m-Lf)\<CT fUs-Lhyds 

Jo 

so, thanks to Grownwall's lemma, — ^jUtv = for all t G [0,T]. 
Moment estimate 

There exist a constant Ct depending only on T such that: 

k i\3 



□ 



supE sup ifitA) ] < Ct 



(31) 



This result is similar to Lemma 5.3 Indeed, by dividing each term of Equation (17) by 



k, we can apply the same reasoning thanks to the remark at the end of the proof of the 
lemma. 

Tightness for fixed / 

We first show that the family of laws of [jjt , f) is tight, for each / € C|(Af). For fixed 
t < T, since / is bounded, we have 

P(K/z^/)| >i^) <^Cf supe( sup (/i^l) 

,k 



Using (31 ), we conclude that the sequence of the laws of ((/Xj , f))keN is tight. Then denote 
by j4j the finite variation part of (/i^, /). We have 



A 



A=d,b,c,g 



1 



with 



X 



X''{x)f{x)^l{dx) <Cf{fi,i 



^'(^) [ Ix fi(P- + ^-)) Di^^ ^) d^] ^(d^) ^ (/i, 1) 

X 

[f^kuHx,y)fi{dy)] /.(dy)| <C/(^,1)2 

[V/(x) • g^x, k^i) + l 4, {x, k ^)] /i(dx) < Cf 1) . 



X 



X 
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Hence for any sequence (T„)„gpj of stopping times bounded by T, 



EIA 



At I < Ct 



using (31) again. Similarly, we can prove 



Tk+e 



.k, • -ce 



k 



thanks to (25) and mlv. According to the Aldous-Rebolledo criteria, this ensures the 



tightness of the laws of (/i'^, /). 



Tightness 



The last result allows us to apply Theorem A.l, to conclude that (/i^)^^^ is relatively 
compact in D([0,T], (A4ir(Af), vague topology)). Denote by fi. the limit of any convergent 
subsequence [fi^ )k'm and notice that fi. is a.s. strongly continuous. Indeed, we can easily 
check that 



sup sup 

o<t<Tfecl^x) 

\\f\\oc<l 



{^4'J)-{^4-.f) 



< 



1 



holds true, which in turn, implies (40). 



Caracterization 



We expect the martingale and residual terms to vanish when passing to the limit in (23). 



This would yield the characterization (29) provided some continuity property holds. 



Lemma 6.4 Let jj, taking values in C{[Q,T], Aipii^)) be the limit of any convergent sub- 
sequence (/i^ )fe'eN- Then for any fixed t <T, f ^ C'^{X), the function: 



MO = {Ct,f)-{Co,f)- / if{Cs)ds 

Jo 



is continuous at point onO{[0,T],Mp{X)), for all uj a.s. 



Proof As is a continuous process, from the characterization of the Skorohod met- 
ric (see [m Proposition 6.5 Ch. 3]) we get: (^^ , /) tends to {iJ-t, f) for all t and / 
bounded/continuous. Then: 

|^t(/') - %{^)\ < K^f - f,t, f)\ + K^g' - ^0, /)| + E f \^"f(^'s) - ^""/(/is)| ds 

A=d,b,c,g'^° 

Both "d" and "b" terms are of the form (/, /i^' — /Ug) which converges to 0. Consider the "c" 
term, we have tf{^^g ) — >• I''f{^s) because the function Cs ^ f f u{x, y) f{x) (^s(g)Cs(dx, dy) 
is continuous (for the weak topology). 
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For the "g" term: 



Ix 



+ 



Qf{x,lJ.s) [/if (dx) - /is(dx)] 



X 



with g{x,Q = Vf{x)-g{xX) + 



1 d^f(x) 

2 dxi dxf. 



-a^g,{x,C). Thanks to the Lipschitz continuity 



(27), the first term of the RSH of the previous inequahty is bounded by: 



xex 

which tends to as —^Hs- The first term also tends to 0. 
Note that 



.k' 



.k' 



□ 



where Z^; and are given by (21) and (24) respectively. is a centered martingale, we 



now check that its quadratic variation tends to 0. Indeed, Using (25) and (31), we have 







which shows that Z^ goes to in L^. Similarly, since 



<C/((^,l) + (/X,l)2) , 



A=c,i 



we conclude that goes to in L} thanks to (31) and Lebesgue's Theorem. Finally, for 

CGO([o,r],A^p(A')), 



\^t{C)\<Cf,t sup ((Cs,l) + (Cs,l) 

Q<s<t 



so that {^t{f^^'))k' is uniformly integrable by (31). Lebesgue's theorem and a.s. continuity 
oi ^ at fi yields 

= limE|^f(/')| =E|^'t(;u)|. 

k' 



Using asumption (28), we conclude that /x converges to a process which is a.s. equal to 



the deterministic unique solution of (29). 



Weak convergence 

So far, we have established the convergence of n'^ to ^ in ID)([0, T], {M-p{X), vague topology)). 
The extension to the weak topology is achieved thought [22[ Th. 3]. Indeed, we only need 
to prove that the total mass (/i'^, 1) converges in law to (^., 1) in D([0, T], M), provided that 
the limitting process is continuous. But this is a particular case of the work done with 
/ = !• 



This ends the proof of Theorem 6.1 
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7 Conclusion 



Using the approach described in this paper, most IBMs could be rewritten as measure- 
valued Markov branching/diffusing processes. This allows first an insight into the Monte 
Carlo procedure adapted to the situation, and second to derive a model at macroscopic 
level as a limit model on various asymptotic situations. 

The model presented in this work assumes that there is no limitation for resources. For 
concrete applications, this model should be coupled with models of the resource dynamics. 
It is also important to devote more work to the case of ecosystems in the presence of several 
species. 

Prom the simulation point of view, the direct microscopic IBM approach is limited 
to relatively small population sizes. For realistic scenarios it could be possible to use a 
coupled approach: species of importance with small population sizes, could be modeled 
with IBMs: while less important species with large population sizes, could be modeled 
through coarser macroscopic models. 

Appendix 

A Weak convergence of sequence of measure-valued pro- 
cesses 

Consider the space D([0,T], A1_F(Af)) of cddldg functions with values in the space Aipi'^) 
of finite positive measures. The space I]){[0,T],Mf{'^)) is endowed with the classical 
Skorokhod topology. The space A4f{X) will be equipped with the weak or the vague 
topologies: the weak (rcsp. vague) topology corresponds to the convergence Vn "'""'^> 
(resp. Vn ^'"^"°> I') defined by: 

{un, (t>) > {v, (t>) G Cb{X) (resp. Ck{X)) 

n—^oo 

where is the set of bounded continuous functions from X to TZ (resp. Cxi'^) is 

set of continuous functions from A' to M with compact support). Hence the Skorokhod 
topology on D([0,r], A^i?(A')) could be associated to the weak or to the vague topology 
on Mf{X)- Unless otherwise specified, D([0,T], A^i;'(A')) will be equipped with the weak 
topology. 

Consider a sequence of random variables {fJ,^)keN which take values in D([0, T], A4f{X)) 
(i.e. {fJ-t)te[o,T] is a cddldg process which takes values in A4f{X))- Here we want to sum- 
marize the conditions that ensure that /j,^ converges in law, that is converges weakly 
where Q'^ denotes the law of on I]){[0,T],M.f{'^))- If Q is the limit law and /i is a 
process with law Q, then 

Q'^ ^ Q (32) 

A;— >-oo 
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means E^{fL^) — )• E$(/x) for all function $ : D([0, T], A^ir(Af)) i-)- M continuous and 
bounded. Hence ( p2| ) is: 



law 

^/^ 

fe— >oo 

This property is sometimes called weak convergence of the processes fx^ toward [2T] . 

The study of the convergence of is usually accomplished in two steps [21] : 

(si) The tightness step: First one shows that the sequence is relatively compact]^ It is 
equivalent to the uniform tightness property: for every e > there exists a compact 
set ICe in D([0, T],Mf{X)) such that inf^ Q''{}C^) > 1 - e (Q. By extension, we will 
say that the sequence iJ,^ is tight if the sequence of their laws is tight. 

(S2) The limit uniqueness step: One proves that there is a set of properties that are 
fulfilled by the limit of any convergent subsequence of Q^. By showing that 
only one law Q satisfies this set of properties, we both prove that the sequence Q'^ 
converges and characterize the limit law Q with the associated limit process /i. This 
set of properties is usually expressed in terms of a martingale problem. 

The tightness step (si) 

The tightness step (si) is achieved through the: 

Theorem A.l ( |14L Theorem 9.1 Ch. 3]) Suppose the the following compact contain- 
ment condition holds: for every e > there exists a compact set of Aipi^) such that 

^f^^fJ't Ke forO<t<T^ >l-e. (33) 

Let Q be a dense subset of CbiMpi^)',^) for the topology of uniform convergence on 
compact sets. Then {lJ.'^)keN is relatively compact if an only if (/(/u'^))fegN is relatively 
compact as a family of processes in D([0, T],M) for any f €z Q. 

An example of such a set is given in [25, Th. 2.1] by G = {/^ ; (j) G &}, with 
f(t>{i^) = W, where is any set dense in Cb{X). 

For any given / G 0, to check that: 

is tight in D([0, T],M). To this end, when (^'^)tg[o,r] is a semimartingale, we make use of 
the Aldous-Rebolledo criteria that can be found in [211 Cor. 2.3.3]: 

■^i.e. with compact closure, it is equivalent to the fact that from any subsequence one can extract a 
convergent subsequence. 

^According to Prohorov's theorem: Tightness implies uniform relative compactness, and it is equivalent 
because J]){[0,T], Mf{X)) with the Skorohod metric is complete and separable (see [141 Theorem 2.2 Ch. 
3]). 
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Theorem A. 2 (Aldous-Rebolledo criteria) Let {Y^^)k£N be a sequence of real valued 
semimartingales with cadlag paths. If: 

(i) For each fixed t G [0, T], (l^^)^^^ is tight. 

(a) Let = + where is a finite variation process and is a locally square- 
integrable martingale. Suppose that for any given sequence of stopping times t^, 
hounded by T, for each e > there exists 6 > and k^: 



sup sup P _ > e < £ , (34) 

fc>feoee[o,<5] ^ ^ 

sup sup ¥(\{M^)r,+e - {M^)rA >e)<e. (35) 
k>ko ee[o,<5] ^ ^ 

Then the sequence {Y^)k^^ is tight. 

In order to get a compact set Kg, of A4.p{X) for the compact containment condition 



(33), we can consider a set of the form: 

K,'^[^,&MF{x)■A^^^)<Ce}. (36) 

When X is not compact, Kf, is compact for the vague topology [TTl Th. 13.4.2] and not 
compact for the weak topology. When X is compact, the vague topology and the weak 
topology are the same, and Ke is compact. 

(i) A first approach proposed in [13] is to consider the compactification X X. Then 



Kg, the set ( 36 ) where X is replaced hy X, is weakly compact and we get the tightness 
on D([0,T], (A^_f(A'), weak topology)). To conclude, we need to prove that the limit 
process lies on A^i?(Af), i.e. Q is supported by O{[0,T], M.p{X)). 

(a) A second approach proposed in [25j is to work with the vague topology: we first prove 
tighness on D([0, T], (A^i?(^Y), vague topology)) and a the associated convergence 
result. To conclude, according to [22, Th. 3], if: 

Jaw^ ^ ID)([0,r], {Mf{X), vague topology)) , (37) 

k—^oo 

J!^(^ ,l) in D([0, T];M), (38) 
^. is a process on C{[0,T]; M.f{X)) , (39) 



then (32) holds true. 



According to [2^ Prop. 3.26], to ensure ( |39[ ) we prove that for all e > 0: 

lim P ( sup |(^,^ /) - (m^ , /) I > = (40) 

for all / G C^,{X). 
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The limit uniqueness step (S2) 

For the step (32): let Q be the Hmit of a any convergent subsequence Q^' and /i. a process 
with law Q. If we prove that /i. satisfies an integro-differential equation that admits a 
unique solution, then ^ Q. 

Usually, Q is characterized as a solution of a martingale problem which appears to be 
the limit of martingale problems satisfied by Q^. At this level, an argument allowing us 
to pass to the limit in the martingale problem should be invoked as the one proposed in 
[14:\ Ch. 4 Th. 8.10]. Then the uniqueness of the solution of the limit martingale problem 
could be obtained by a duality argument as proposed in [131 Section 1.6] or [Tl', p. 188]. 
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